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ABSTRACT 

Using conditional luminosity functions (CLFs) which encode the luminosity distri- 
bution of galaxies as a function of halo mass, we construct a halo model of IRAS 
galaxies selected at 60 /im. An abundance matching technique is used to link galaxy 
luminosity to the host halo mass. The shape of the mass - light relation at 60 /j,m is 
different from those derived at r-, K- and B-band. This is because the 60 /im LF can 
not be fitted by a Schechter function with a sharp exponential cutoff above L*. We 
then seek the parameters in the CLFs that best fit the 60 /im LF and power spectrum. 
We find that the predicted galaxy bias as a function of Lqq from the best-fit model 
agrees well with the clustering measurement in different luminosity bins. At the faint 
end of the LF where quiescent star-forming galaxies dominate, most IRAS galaxies 
are central galaxies in halos of M > 10 10 h^ 1 Mq but a non- negligible fraction are 
satellites typically hosted in more massive halos. The majority of IRAS galaxies with 
^60 ^ 10 10 h~ 2 Lq are M82 type starbursts which are central galaxies hosted in halos 
of M > 10 12 ' 5 h~ l Mq. In comparison, optical galaxies generally reside in much more 
massive halos. The rate of change in Lqq (an indicator of recent star formation) as 
a function of halo mass at M > 10 125 h~ 1 M & is much larger that dL op ti ca \/dM or 
g^nir/ 'dM indicating the existence of physical mechanisms which are very efficient 
in converting cold gas into stars, possibly dynamical effects arising from interactions 
or mergers. We further calculate the space density of major mergers for halos mas- 
sive enough to host ultraluminous infrared galaxies (ULIRGs) using the mean merger 
rate derived from the Millennium simulations. Compared to the space density of local 
ULIRGs, it implies that either the majority of major mergers at z ~ do not lead to 
ULIRGs or the ULIRG phase is relatively short. 

Key words: infrared: galaxies - galaxies: haloes - galaxies: luminosity function, mass 
function - galaxies: starburst - galaxies: spiral - large-scale structure of Universe. 



1 INTRODUCTION 

Understanding the relative distribution between different 
populations of galaxies and the underlying dark matter is 
one of the most important goals in modern astrophysics and 
cosmology. In the local Universe, large galaxy spectroscopic 
surveys covering thousands of square degrees (e.g. SDSS 
and 2dFGRS) reveal how galaxies are distributed spatially 
and how that varies with intrinsic galaxy properties such as 
colour, luminosity and spectral type. A lot of progress have 
also been made in mapping galaxy distribution in the distant 
Universe albeit over much smaller areas (e.g. COSMOS and 
DEEP2). On the other hand, numerical simulations give us 
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a fairly good understanding of the abundance and cluster- 
ing of dark matter haloes. There are in general three differ- 
ent approaches to connect galaxies with dark matter halos. 
Hydro-dynamical simulations attempt to predict galaxy for- 
mation and evolution from first principles but it is computa- 
tionally too expensive to carry out cosmological simulations 
with enough resolution and volume (e.g. Katz et al. 1992; 
Evrard et al. 1994; Frenk et al. 1996; Weinberg et al. 1997; 
Springel & Hernquist 2003; Springel 2005; for a review see 
Springel 2010). In the foreseeable future, it seems infeasible 
to use hydro-dynamical simulations to explore a reasonable 
range of parameter space related to baryonic physics. Semi- 
analytic models (SAMs) also attempt to simulate the for- 
mation and evolution of galaxies in an a priori fashion but 
with highly simplified analytic prescriptions to approximate 
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star formation and feedback processes (e.g. Kauffmann et 
al. 1999; Somerville & Primack 1999; Benson et al. 2000, 
2001; Springel et al. 2005; Somerville et al. 2008). It is com- 
putationally inexpensive which means that it is possible to 
use MCMC techniques to find the best-fitting parameters 
to reproduce observables such as galaxy luminosity func- 
tion (LF) , mass function and clustering. However, there are 
many degeneracies between various parameters in SAMs. 
More importantly, many physical processes (such as differ- 
ent feedback mechanisms) are still poorly understood which 
results in large uncertainties in various recipes used to follow 
galaxy formation and evolution. In recent years, an empir- 
ical approach based on virialized halos has been developed 
to circumvent these physical uncertainties. The basic halo 
model of the large-scale distribution of galaxies is a statisti- 
cal approach that links galaxies with their host dark matter 
halos (e.g. Jing et al. 1998, 2002; Peacock & Smith 2000; 
Seljak 2000; Scoccimarro et al. 2001; Berlind & Weinberg 
2002; Cooray & Sheth 2002; Zheng 2004; Zehavi et al. 2004, 
2005; Collister & Lahav 2005; Tinker et al. 2005; Skibba 
et al. 2007; Brown et al. 2008). The halo occupation distri- 
bution (HOD) gives the the number of galaxies of certain 
intrinsic properties (such as luminosity, colour) in a halo of 
mass M. Combined with an assumption about the spatial 
distribution of these galaxies in the halo, the abundance and 
clustering of haloes can be used to derive those of galaxies. 
One of the disadvantages of the halo model is that it does 
not directly provide a physical understanding of the result- 
ing halo occupation function. 

In this paper, we apply an important extension of the 
basic halo model, the so-called conditional luminosity func- 
tion (CLF), to local IRAS galaxies selected at 60 jj,m. CLF 
<E>(L|M) describes the number of galaxies as a function of 
luminosity and halo mass (Yang et al. 2003, 2005; Cooray 
& Milosavljevic 2005; Cooray 2006). Using CLF, one can 
naturally address the clustering and abundance of galaxies 
as a function of luminosity which is difficult to study us- 
ing the basic halo model. Both the basic halo model and the 
CLF have been successfully applied to local galaxies selected 
from large spectroscopic surveys such as SDSS and 2dFGRS 
and distant galaxies selected from DEEP2, COMBO-17 (e.g. 
Zehavi et al. 2005; Cooray & Milosavljevic 2005; Cooray 
2006; Phleps et al. 2006; Coil et al. 2006; Zheng et al. 2007). 
For example, the halo model can successfully fit the statis- 
tically significant deviations from the canonical power-law 
behaviour in galaxy correlation functions which in turn con- 
strains parameters related to galaxy formation and evolu- 
tion processes (e.g. the halo mass scale for hosting a central 
galaxy and satellites). However, so far there are no avail- 
able halo models of local IRAS galaxies. Infrared luminosity 
is a good indicator of star formation activity with is not 
affected by dust obscuration. Understanding how infrared 
galaxies populate dark matter haloes could give us impor- 
tant insights on the relation between star formation and 
environment and the possible mechanisms for triggering / 
quenching star formation which may preferentially operate 
on different mass scales (e.g. ram pressure stripping, major 
merger). 

This paper is organised as follows. In Section 2, we 
briefly summarise the ingredients in the CLF halo model. 
Further details can be found in Cooray & Sheth 2002, Yang 
et al. 2003, 2005, Cooray 2006. In Section 3, firstly we de- 



rive the central galaxy luminosity - halo mass relation and 
the total galaxy luminosity - halo mass relation using an 
abundance matching technique. We then constrain the two 
free parameters left in the CLF model using the 60 /jm LF 
and the power spectrum of IRAS galaxies selected from the 
Point Source Catalogue Redshift (PSCz) catalogue (Saun- 
ders et al. 2000) . In Section 4, the best-fit CLF model is used 
to predict the effective galaxy bias as a function of 60 /xm lu- 
minosity which agrees well with the clustering measurements 
of IRAS galaxies in different luminosity bins. In Section 5, 
we derive the comoving number density of major mergers 
which is then compared to the measured number density of 
ultraluminous infrared galaxies (ULIRGs). Discussions and 
conclusions are presented in Section 6. Throughout the pa- 
per, we use a spatially flat ACM cosmology with 57a = 0.7, 
Q m = 0.3, the slope of the primordial power spectrum 
n s = 1 and the overall normalisation of the power spec- 
trum erg = 0.8. In addition, we build our model at redshift 
z = and redshift evolution is not considered. 



2 CONDITIONAL LUMINOSITY FUNCTIONS 
(CLF) 

The CLF gives the average number of galaxies that resides 
in a halo of mass M as a function of luminosity. The CLF is 
divided into two terms corresponding to central and satellite 



$(L|M) = $ ccn (L|M) +$ sat (L|M) 



(1) 



The first term is the number of central galaxies as a function 
of luminosity and halo mass which is assumed to obey a log- 
normal relation, 
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The mean luminosity L C (M) encodes the relation between 
central galaxy luminosity and halo mass and cr cen represents 
the dispersion in this relation. The second term corresponds 
to the luminosity distribution (assumed to be a power-law) 
of satellite galaxies which appear when the total luminosity 
inside a halo exceeds the central galaxy luminosity, 



$ sat (L|M) = A{M)L 1 . 



(3) 



The normalisation $(M) and A(M) are such that 
/ & ccn (L\M)LdL = L C (M) and / $ sat (LjM)WL = 
Ltot(M) — L C (M). The galaxy LF can then be obtained by 
integrating CLFs over the halo mass function, 



*(L) = J 



dMn(M)$(L\M). 



(4) 



Here n(M) represents the mass function of dark matter halos 
and we use the formalism of Sheth & Tormen (1999), 
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where p is the mean matter density of the Universe at z — 0, 
a(M) is the linear rms mass fluctuation on mass scale M, S c 
is the critical overdensity required for collapse at z = 0, v = 
8 c /a(M) , v' = yf{a)v, a = 0.707, q = 0.3 and A w 0.322. 
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Figure 1. The relation between dark matter halo mass and galaxy luminosity. Different solid lines represent the relations between halo 
mass and central galaxy luminosity at different wavebands (black solid line - IRAS 60 /an; red solid lines - SDSS r-band; green solid line 
- rest B-band; blue solid line - rest K-band). Different dashed lines represent the relation between halo mass and total galaxy luminosity 
( black dashed line - IRAS 60 ^m; blue dashed line - rest K-band. The horizontal solid lines mark the locations above which the 60 fim 
LF is dominated M82/Arp220 type starburst. The horizontal dot-dashed lines indicate L, at different wavebands. 



Similarly, the galaxy power spectrum can be decom- 
posed into two clustering terms in the halo model 



P gal (fc|L) = Pi h (fc|L) + P 2h (k\L). 



(7) 



The 1-halo contribution comes from galaxies which reside in 
the same dark matter halo and it can be expressed as 
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+ 2$ con (L|M)$ sat (L|M)« g (fc|M)]dM, 



(8) 



where n(L) is the mean number density of galaxies as a 
function of luminosity, 



n(L) = / dMn(M)Q(L\M) 



(9) 

and u g (fc|M) is normalised Fourier transform of the galaxy 
density distribution within a halo of mass M, 

2 sin kr p g (r\M) 



u g (k\M) = / dr4nr 

'o kr 



M 



(10) 



We assume that the galaxy density distribution within the 
dark matter halos p s (r\M) can be described by an NFW pro- 
file (Navarro, Frenk & White 1997) truncated at the virial 
radius, 

M vir = ^7n^ ir( oA vir , (11) 

where the virial overdensity is A v i r m 337 (Bullock et al. 
2001). The 2-halo contribution comes from galaxies which 



reside in separate dark matter halos, 
fth(fc|L) = P lin (k) 



dMn(M)b(M) 



<b(L\M) 



u g (k\M) 



(12) 



Here P bn (k) is the linear dark matter power spectrum and 
b(M) is the bias factor for dark matter halos of mass M. 
In this paper, we use the bias factor derived by Sheth, Mo 
& Tormen (2001) based on the ellipsoidal collapse model. 
If one is interested in the galaxy power spectrum above 
a certain luminosity threshold or over a wide luminosity 
range, one can simply integrate the CLF $ ccn (L|M) and 
$ sat (L|M) over luminosity to obtain the halo occupation 
numbers N ccn {M) and iV 8at (M) and replace the CLF with 
the halo occupation numbers in Eq. 8, 9 and 12. 



3 OBSERVATIONAL CONSTRAINTS: 

LUMINOSITY FUNCTION AND GALAXY 
POWER SPECTRUM 

Following Vale & Ostriker (2004) , we derive the relation be- 
tween the mass of the dark matter halo/subhalo and the lu- 
minosity of a galaxy hosted in it using an abundance match- 
ing technique. The inputs in this approach are the halo mass 
function n(M) (Eq. 5 and 6), the mass distribution of sub- 
halos in a halo of a given mass and the galaxy LF. The mass 
distribution of subhalos in a parent halo of mass M is 
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N(m\M)dm = A 



x/3M 



exp 



dm 



x/3M ) x/3M ' 



(13) 



where a = 1.91, j3 = 0.39, 7 = 0.18, x = 3 and 4 = 
7/[/3r(2-a)] (De Lucia et al. 2004; Weller, Ostriker & Bode 
2005). The total mass in these subhalos is a fraction of the 
parent halo mass, x-yM. The global subhalo mass function 
can thus be obtained by integrating over the halo mass func- 
tion.. 



n sh (m) 



N(m\M)n(M)dM. 



(14) 



The halo mass - central galaxy luminosity relation can be 
derived by matching the number density of galaxies above 
a certain luminosity and the number density of haloes / 
subhaloes above a certain mass threshold, 

&(L)dL= [n{M) +n sh (M)]dM. (15) 

lL JM 

We use the present-epoch 60 fim LF from Saunders et al. 
(1990) which is a Gaussian combined with a power law, 
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L* = 10 s ' 47 hr 1 L @ . Note that f>(L) is the LF per decade 
in luminosity. There are two implicit assumptions in this 
approach: (1) the luminosity of a galaxy hosted in a halo 
is an increasing function of the halo mass; (2) Each halo or 
subhalo can only host up to one galaxy in total. The halo 
mass - total luminosity relation can be expressed as follows 



Ltot(M) = L(M) + L(m)N(m\M)dm, 



(17) 



which means the total luminosity in a halo of mass M is the 
sum of its central galaxy luminosity and satellite luminosi- 
ties in all its subhalos. 

Fig. [T] shows the derived mass - luminosity relation for 
IRAS galaxies selected at 60 /im. In halos less massive than 
1O 12 /i _1 M , the mass - total luminosity relation is essen- 
tially the same as the mass - central galaxy luminosity re- 
lation. Compared with mass - light relations at other wave- 
bands, e.g. SDSS r-band, rest B-band and rest K-band, the 
60 ;Um mass - light relation seems to be distinct in the sense 
that it can not be fitted by the general fitting formula, 



L C (M) = L - 



(M/Mi 



'[&+ (M/Mi)**} 1 /*' ^ 

This is because the 60 ^m LF does not show a sharp ex- 
ponential cutoff above L* and therefore is not well fitted 
by the Schechter function = $»(L/L»)~ Q exp(-L/L») 

(Schechter 1976) which accurately describes the LFs at r-, 
B- and K-band (e.g. Blanton et al. 2001, 2003; Norberg et 
al. 2002; Smith et al. 2009). In other words, there are rela- 
tively more luminous infrared galaxies than luminous optical 
galaxies. 

The knee of the LF L* is shown in Fig. [1] for different 
wavebands. It is immediately obvious that galaxies with lu- 
minosities around L* at r- and K-band reside in halos of 
similar mass while L, galaxies selected at 60 fim are hosted 
by much less massive halos (roughly by over an order of 
magnitude). This is consistent with the relative bias be- 
tween these different tracers of the large-scale structure de- 
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Figure 2. The IRAS PSCz galaxy power spectrum (empty dia- 
monds) measured in Hamilton & Tegmark (2002) compared with 
the predicted power spectrum from the best-fit CLF model. 



rived from clustering measurements (e.g. Peacock & Dodds 
1994; Fisher et al. 1994; Mann et al. 1996; Hawkins et al. 
2001; Zehavi et al. 2002, 2005; Ma et al. 2009). In Fig. [[] 
we also show the luminosity scales Lao = 10 10 K~ 2 Lq and 
Lbo = 10 11 ' 4 Ii~ 2 Lq above which the 60 /im LF is domi- 
nated by M82 type (interaction-induced) and Arp220 type 
(merger-induced) starbursts, respectively (Wang & Rowan- 
Robinson 2010). At luminosities below L 60 = 10 10 /t~ 2 Le, 
the 60 /jm LF is dominated by cirrus type quiescent star- 
forming galaxies. In the luminosity range where M82 and 
Arp220 type starbursts dominate, the rate of change in Leo 
as a function of halo mass M is much larger than the rate 
of change in Z/ op ti C ai or Lnir. Since the 60 /im luminosity 
can be used as a proxy for star formation rate, it implies 
that in halos with masses > 1O 12 5 /i _1 M , there are some 
physical mechanisms which are very efficient in converting 
cold gas into stars. From numerical simulations, groups have 
a virial mass of ~ 10 13 Mq, while clusters have virial 
masses of 10 14 (poor) - 10 15 h,- 1 Mq (rich). So, one possibil- 
ity is that dynamical effects such as gravitational interaction 
and merging activity in these more massive halos trigger a 
frantic phase of star formation. Observationally, most local 
luminous infrared galaxies show signs of interaction or merg- 
ing activities at various stages (e.g. Sanders & Mirabel 1996; 
Clements et al. 1996; Murphy et al. 1996; Farrah et al. 2001; 
Veilleux, Kim & Sanders 2002). Numerical simulations have 
also shown that tidal forces in interacting / merging galax- 
ies are very efficient at creating centrally concentrated gas - 
fuel for nuclear starburst and active galactic nuclei activity 
(e.g. Mihos & Hernquist 1994a,b, 1996; Springel 2000; Cox 
et al. 2006; Di Matteo et al. 2007). We will return to this 
point in Section 5. 

Now we have derived the relation between central / to- 
tal galaxy luminosity and halo mass, the only free parame- 
ters left in CLF are the power-law slope 7 in the luminosity 
distribution of satellite galaxies and the dispersion a con in 
the luminosity distribution of central galaxies. We vary both 
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4 TESTING THE BEST-FIT MODEL 
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Figure 3. The halo occupation numbers as a function of halo 
mass. The solid lines correspond to satellite galaxies and the 
dashed lines correspond to central galaxies. The red, green and 
blue curves represent the halo occupation numbers for galaxies in 
the luminosity range [10 s , 10 9 ], [10 9 , 10 10 ] and [10 10 , 10 11 ] h~ 2 
Lq, respectively. 



<T ccn and 7 over a wide range and search for the best-fit val- 
ues that minimise \ 2 = X 2 (®( L )) + X 2 (P( k ))- For the LF 
we use the non-parametric maximum likelihood esti- 
mates from Wang & Rowan-Robinson (2010). For the galaxy 
power spectrum P(k), we use the real-space power spectrum 
of the IRAS PSCz survey from Hamilton & Tegmark (2002) . 
We find that the best-fit parameters are 7 = —2.96 and 
(icon = 0.022. In Fig. we compare the IRAS PSCz power 
spectrum with the predicted power spectrum from the best- 
fit model. The 1-halo clustering term is sensitive to 7 - the 
power-law slope in the satellite luminosity distribution. The 
transition from the 2-halo term to the 1-halo term occurs 
at k ~ 1 h Mpc~ x . Fig. [3] shows the halo occupation num- 
bers for central and satellite galaxies as a function of halo 
mass. At a give halo mass scale, the number of satellites 
decreases rapidly as Leo increases. At a given luminosity 
scale, the number of satellites increases rapidly as a func- 
tion of halo mass. It indicates that a significant fraction of 
low-luminosity IRAS galaxies are satellites hosted in massive 
dark matter halos. In Fig. |3J the predicted LF and the con- 
tribution from central and satellite galaxies from the best-fit 
CLF model are compared to the measured LF. The 60 )im 
LF is dominated by central galaxies at all luminosities apart 
from the faintest luminosity bins (Leo ~ 10 7 ' 5 h~ 2 Lq), sim- 
ilar to the finding at other wavebands (Cooray et al. 2005). 
The shape of the LF at the bright end is most sensitive to 
(j ccn . At luminosities below Leo = 10 10 h~ 2 Lq where the 
60 fim LF is dominated by cirrus type galaxies, the frac- 
tion of satellite galaxies increases rapidly as Leo decreases. 
This implies that many satellites in massive halos could have 
on-going quiescent star formation. 



Using the real-space galaxy auto-correlation functions, 
Wang & Rowan-Robinson (2010) measured the bias factor 
as a function of 60 /jm luminosity using galaxies from the 
Imperial IRAS-FSC Redshift Catalogue (IIFSCz; Wang & 
Rowan-Robinson 2009). In this section, we test our best- 
fit CLF model using the luminosity dependence of galaxy 
clustering. From the best-fit CLF model, we can derive the 
effective bias factor as a function of luminosity, 

$(L|M) 



b e (L) = / dMn(M)b(M)- 



(19) 



In Fig. [5l we plot the predicted bias factor as a function of 
60 /im luminosity. A good agreement (x 2 e d = 1-1) is found 
between the model prediction and the measured bias fac- 
tors. The dashed line in Fig. [5] shows galaxy bias as func- 
tion of luminosity for optical galaxies derived from mea- 
surements of the SDSS power spectrum (Tegmark et al. 
2004), b g {L)/b g (L„) = 0.85 + 0.15L/L, - 0.04(M — M»). 
We use L* = 10 8 47 h~ 2 Lq for IRAS galaxies and L* = 
10 10 h~ 2 Lq for optical galaxies. Compared with the bias de- 
rived from optical galaxies, the clustering properties of IRAS 
galaxies show a rather mild luminosity dependence around 
L, . This is because IRAS L» galaxies reside in low mass ha- 
los where the linear bias factor does not change rapidly as a 
function of halo mass (e.g. Sheth & Tormen 1999; Sheth, Mo 
& Tormen 2001; Seljak & Warren 2004). The rise of galaxy 
bias towards lower Leo is caused by an increasing fraction 
of low luminosity IRAS galaxies being hosted as satellites in 
more massive halos. In addition, we do not find any evidence 
for merger bias which is a tendency of recently merged sys- 
tems to be more strongly clustered on large scales than typi- 
cal systems of similar mass (e.g. Furlanetto & Kamionkowski 
2006; Chapman et al. 2009). 



5 ULIRGS AND MAJOR MERGERS 

Local ultraluminous infrared galaxies (ULIRGs) with in- 
frared luminosities Lm ^ 10 12 h~ 2 Lq are believed to be 
gas-rich major merger systems which are evolutionally con- 
nected with QSOs (e.g. Mihos & Hernquist 1996; Sanders & 
Mirabel 1996; Moorwood 1996; Lonsdale, Farrah & Smith 
2007; Wang et al. 2010). Combining the two Millennium 
simulation data, Fakhouri O. et al. (2010) derived the fol- 
lowing fitting formula for the dimensionless mean merger 
rate (i.e. the mean number of mergers per halo), 



dN m 
d£dz 
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10 12 M, 



£ exp 
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where M is the descendant halo mass, £ is the progeni- 
tor mass ratio (i.e. M1/M2 with Mi < A/2), (a, /J, 7, 77) = 
(0.133,-1.995,0.263,0.0993) and = (0.0104,9.72 x 



10" 



The total space density of z w major mergers n n 



for halos hosting ULIRGs can be calculated as 



1 dN 
1/3 d£ dz 



[M,S,z = Q)dtW, (21) 



where Mulirg, the average halo mass scale hosting infrared 
galaxies with Lir = 10 12 h~ 2 Lq, is derived from the mass - 
light relation in Fig.[T]using an empirical ratio of Lir/L60 = 
2.2T0.2. Using a volume- limited sample of ULIRGs from the 
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Figure 4. The symbols are the 60 fim luminosity function of IRAS galaxies derived from a non-parametric maximum likelihood method 
(Wang & Rowan- Robinson 2010). The red solid line corresponds to the contribution from central galaxies while the blue solid line 
corresponds to the contribution from satellites. The black solid lines shows the sum of the two. The red dotted and dashed lines shows 
the contribution of central galaxies in different halo mass ranges indicated by the numbers shown under each curve. The amplitude of 
the red dotted and dashed lines have been reduced by a factor of 100 to avoid overcrowding. 



IIFSCz at z 0.1, we derive the co-moving space density of 
local ULIRGs to be around nuLiac = 2.2 x 10~ 7 h 3 Mpc" 3 . 
Compared to the theoretical estimate in the same redshift 
range n mm = 1.1 x 10 -6 h 3 Mpc -3 , it implies that the ma- 
jority of major mergers are not observed as ULIRGtfl This 
could be due to either a large fraction of the major merger 
events at z ~ are dry mergers rather than wet mergers or 
the ultraluminous infrared phase is relatively short-lived so 
that only a small fraction are observed. There is still a lot of 
debate on the importance of dry mergers in the formation 
of massive galaxies. Some numerical studies show that dry 
mergers dominate the overall merger rate at z ~ (e.g. Naab 
et al. 2006; Khochfar & Burkert 2003, 2005; Khochfar & Silk 
2009) . The short duration of the ULIRGs phase is probably 

1 There are a few caveats here. Firstly, a halo - halo major merger 
may not be the same as a galaxy - galaxy major merger which 
is believed to trigger the ULIRG phase. However, it is not clear 
which mass ratio definition is the correct one to use, stellar mass, 
baryonic mass, dynamical mass or halo mass. Secondly, the con- 
version from the halo - halo merger rate to the galaxy - galaxy 
merger rate is highly uncertain mainly due to the uncertainties 
in simulating baryonic physics. For example, by adopting differ- 
ent merger timescales (i.e. the time delay between a halo - halo 
merger and the subsequent galaxy - galaxy merger) , the predicted 
galaxy - galaxy merger rate can differ by a factor of ~ 2. Hopkins 
et al. (2010) investigated various error sources which in combi- 
nation can lead to order-of-magnitude variation in the predicted 
galaxy-galaxy merger rates. 



more likely to be the cause of the difference between tiulirg 
and n mm from a variety of arguments. Stellar population 
synthesis models show that the UV and optical spectra of 
ULIRGs are consistent with star formation bursts of ~ 10 7 - 
10 8 years old (e.g. Canalizo & Stockton 2000a,b, 2001; Far- 
rah et al. 2005; Rodriguez Zaurin et al. 2008, 2009). Similar 
numbers come from gas depletion time estimates. Assuming 
the bulk of infrared luminosity arises from star formation, 
the gas depletion time is M gas /SFR ~ 10 7 — 10 8 years (e.g. 
Sanders et al. 1988; Solomon et al. 1997; Downes & Solomon 
1998; Evans et al. 2002; Papadopoulos et al. 2008). 



6 DISCUSSION AND CONCLUSION 

We present a halo model of local IRAS galaxies selected at 60 
/im using conditional luminosity functions (CLFs) 4>(L|M). 
We match the number of galaxies above a certain luminosity 
threshold to the number of halos and subhalos above a cer- 
tain mass threshold to derive the relation between halo mass 
and the 60 /im luminosity. The implicit assumptions in this 
approach are: (1) The mass - light relation is monotonic; (2) 
every halo / subhalo can only host up to one galaxy. The 
resultant mass - light relation at 60 fj,m is very different from 
those derived at r, K and B-band at the bright end because 
the 60 /an luminosity function (LF) is not well fitted by 
a Schechter function with a sharp exponential cutoff above 
L„. The typical halo mass scales hosting L* galaxies at 60 
/im, r- and K-band are consistent with the measured relative 
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Figure 5. The measured effective bias of IRAS galaxies as a 
function of 60 (im luminosity compared with the predicted lu- 
minosity dependence from the best-fit CLF halo model. The 
dashed line represents the relative bias for optical galaxies de- 
rived from SDSS power spectrum measurements, b g (L) /b g (L t ) = 
0.85 + 0.15L/L, - 0.04(M - M») (Tegmark et al. 2004). We use 
L* = 10 8 - 47 h- 2 L for IRAS galaxies and L» = 10 10 h~ 2 L 
for optical galaxies. 



bias between infrared galaxies, optical galaxies and near- 
infrared selected galaxies. In addition, the rate of change in 
Leo as a function of halo mass is much larger that the rate of 
change in L op ticai or Lnir, indicating the existence of physi- 
cal mechanisms in these halos which are efficient in inducing 
star formation activities. Gravitational interactions and ma- 
jor mergers are likely to be the cause of high star formation 
rates based on morphological studies of LIRGs and ULIRGs. 

Building on the mass -light relation at 60 /xm, we seek 
the parameters in the CLF model that best fit the 60 /jm 
LF and the power spectrum of IRAS PSCz galaxies. In the 
best-fit model, we find that the majority of IRAS galaxies 
with L 60 > 10 10 h~ 2 L Q are M82 type starbursts which are 
central galaxies hosted in halos of mass > 10 12 ' 5 hT 1 Mq. 
The resultant best-fitting model is also tested using the lu- 
minosity dependence of galaxy clustering. A good agreement 
is found between the predicted bias as a function of 1/60 and 
the clustering measurement of IRAS galaxies in different 
luminosity bins. Compared to optical galaxies, IRAS galax- 
ies show a mild luminosity-dependent clustering around the 
characteristic luminosity L* . This is because IRAS L» galax- 
ies are hosted in low mass halos where the halo bias does 
not change rapidly with mass. In addition, we do not find 
evidence for merger bias. 

Lastly, we compare the predicted major merger number 
density with the measured number density of ultraluminous 
infrared galaxies (ULIRGs) . We find that the former is about 
five times larger than the latter. This could mean that the 
majority of the major merger events do not lead to ULIRGs. 
However, this is probably more likely to be explained by the 
relatively short duration of the ULIRG phase. 

In a future paper, we will use CLF halo model to fully 



model the cross-correlation functions between different types 
of infrared galaxies and the cross-correlation functions be- 
tween infrared galaxies from IRAS and optical galaxies from 
SDSS. This work can also be extended to model the relative 
distribution between infrared galaxies detected by Herschel 
and other classes of galaxies in the distance Universe (Am- 
blard & Cooray 2007). It would also be interesting to com- 
pare the CLF model constrained by observational data with 
predictions from semi-analytic modelling. 
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